Critical point shift in a fluid conflned between opposing walls 
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The properties of a fluid, or Ising magnet, confined in a L x oo geometry with opposing surface 
fields at the walls is studied by density matrix renormalization techniques. In particular we focus on 
the effect of gravity on the system, which is modeled by a bulk field whose strength varies linearly 
with the distance from the walls. It is well known that in the absence of gravity phase coexistence 
is restricted to temperatures below the wetting temperature. We find that gravity restores phase 
coexistence up to the bulk critical temperature, in agreement with previous mean field results. A 
detailed study of the scaling to the critical point, as L — > oo, is performed. The temperature shift 
scales as 1/ L"^ , while the gravitational constant scales as 1 /L^^^^ , with yr and j/h the bulk thermal 
and magnetic exponents respectively. For weak surface fields and L not too large, we also observe a 
regime where the gravitational constant scales as i/L'^+yH-^ivr ^/\^_^ jg ^j^g gyj-face gap exponent) 
with a crossover, for sufficiently large L, to a scaling of type . 

PACS numbers: 05.50.+q, 05.70.Fh, 68.35.Rh, 75.10.Hk 



I. INTRODUCTION 

The Ising model has played a central role in the the- 
ory of critical phenomena. It is simple enough that it 
can be studied in great detail (although exact solutions 
are restricted only to a few cases in two dimensions), 
and it can be used to model many interesting physical 
situations. The two phases with opposite magnetization, 
which coexist in an infinite system for temperatures lower 
than the bulk critical temperature Tc, can be thought as 
the two coexisting phases (liquid and vapor) of a simple 
fluid. Besides bulk critical phenomena, which occur in a 
system infinitely extended in all directions, much interest 
has been focused on the study of the critical behavior of 
confined systems, which are of finite extensions in one or 
more directions. 

A few years ago Parry and Evans using a mean 

field Ginzburg-Landau approach, investigated the phase 
diagram of a d-dimensional Ising model in a L x cxd''^^ 
geometry, i.e. confined between two walls separated by a 
finite distance L. They considered opposing walls, where 
one wall favors the "liquid" and the other the "vapor" 
phase; in the Ising language this corresponds to introduc- 
ing surface magnetic fields hi and /12 with opposite sign 
(/iift.2 < 0). It was found 0] that two phase coexistence 
is restricted to temperatures below the wetting temper- 
ature Tu,; Brochard and de Gcnnes had come to a sim- 
ilar conclusion as well These results were later con- 
firmed by extensive Monte Carlo simulations, both in two 
and three dimensions The two dimensional Ising 

model with opposing surface fields was recently solved ex- 
actly by Maciolek and Stecki using transfer matrices 
and pfaffian techniques. The surprising characteristic of 
Parry and Evans' results is that one does not seem 
to recover information about the bulk critical point when 



the limit i — > 00 is taken: for all values of L two phase 
coexistence is restricted to T < T^,. 

In trying to clarify the remarkable properties of this 
system, Rogiers and Indekeu Q| introduced an extra bulk 
field which varies linearly with the distance from the walls 
and which plays the role of gravity in the fiuid. They 
found, in a mean field analysis, that the competing ef- 
fect of opposing surface fields and gravity restores phase 
coexistence up to the bulk critical temperature. 

In the present paper we study the model of Rogiers and 
Indekeu at the lower critical dimension (d — 2) and be- 
yond the mean field approximation to test the validity of 
their conclusion when thermal fluctuations are properly 
taken into account. Our results are based on a density 
matrix renormalization group (DMRG) calculation and 
essentially conflrm the mean held scenario. We also ana- 
lyze the exponents describing the critical point shift along 
the thermal and gravitational held direction and compare 
our flndings with some predictions due to scaling analysis 

"Cli. 



A brief report of this work has been presented in Ref. 
p2| ; here we give a full account of our work and present 
a series of results concerning in particular the analysis 
of the finite size scaling of the gravitational constant for 
low surface fields, where interesting crossover phenomena 
take place. This paper is organized as follows. In Sec.|| 
we introduce the model and in Sec, III briefiy review the 
DMRG technique. In Sec.IV we compare DMRG results 
with exact ones in the absence of a gravitational field. 
In Sec.0 we present the phase diagram as derived from 
DMRG calculations and discuss its main features. In 
Sec.0 we show some magnetization profiles and com- 
pare the DMRG resu lts w ith those obtained by capillary 
wave theory. In Sec. VII we present some conclusions. 
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II. MODEL 

We consider a ferromagnetic Ising model in a L x cxd 
strip with the foUowing Hamiltonian: 



-/ii^SLj+.g^^(2i-l-L)5,,j, (1) 

where J > 0, s^j = ±1, — oo < j < oo and 1 < i < L. 
For simphcity we restrict ourselves to the case of surface 
fields of equal magnitude \hi\. The last term of the right- 
hand side of Eq. (|l]) is a bulk field varying linearly with 
the distance from the walls, which, for convenience, is 
chosen antisymmetric with respect to the center of the 
strip. This field models the effect of gravity on a fluid 
and g plays the role of the gravitational constant. 

In the limit 5 — > one recovers the models studied by 
Parry and Evans jT] ; in two dimensions the wetting tem- 
perature is known exactly and satisfies the equation 



e° [cosh(a) — cosh(6)] = sinh(a), 



(2) 



where a = 2J/T^, and b = 2hi/T^. For /ii = the 
wetting temperature is equal to the bulk critical tem- 
perature Tc = 2J/ln(l + V2)- It decreases monotonically 
with hi and vanishes for hi = J . 

In the rest of the paper we will be primarily interested 
in the competing effect of surface fields and gravity; it is 
clear from the Hamiltonian (|^) that these occur when hi 
and g have the same sign. 



III. DMRG METHOD 

The density matrix renormalization group was intro- 
duced by White Qjl^] to study the ground state proper- 
ties of quantum spin chains. The method is very accurate 
and it has been successfully applied to many one dimen- 
sional quantum problems (for a review see Ref. [p^). 

Exploiting the relation between a d-dimensional quan- 
tum system and & d + 1-classical system, Nishino 
was able to extend the DMRG to two-dimensional clas- 
sical systems. In this case White's algorithm is applied 
to the construction of effective transfer matrices of large 
systems. 

In a transfer matrix approach the partition function 
of a system defined on a L x 00 strip is equal to the 
largest eigenvalue of the so-called transfer matrix Tl @| . 
Numerical calculations are restricted to strips of small 
widths (typically L < 20 for an Ising model) since the di- 
mension of Tl grows exponentially with the strip width. 

In a DMRG calculation one starts from a transfer ma- 
trix of a small system that can be handled exactly. Using 



the information about the thermodynamics of this sys- 
tem one generates an effective transfer matrix of a larger 
system. The strip width grows at each DMRG iteration 
and the spin space is very efficiently truncated to keep the 
dimensionality of the effective transfer matrix controlled. 

Figure l] shows schematically a transfer matrix element 
generated by the DMRG algorithm. The matrix consists 
of block and spin variables indicated in the figure by rect- 
angles and circles respectively. A block, whose states are 
labeled by a variable ^ which can take m possible values 
only, describes approximately a collection of spins. 



c><> 



FIG. 1. Schematic view of a transfer matrix element of a 
strip of width L generated by DMRG. ^ and s label block and 
spin variables respectively, with ^ — l,2...m and s = ±1. 
The total dimension of the matrix is 4m^ x 4m^. The arrow 
denotes the transfer direction. 

Obviously by allowing larger m one obtains more accu- 
rate numerical results; in a typical DMRG calculation the 
accuracy grows very rapidly with m . In the present 
case we found that for strips up to i = 100 a value of 
m = 40 is sufficient to guarantee a very high accuracy of 
the numerical results. In our calculations we have used 
the finite system method described by White in Ref. , 
a version of the DMRG algorithm designed to accurately 
study finite size systems. For more details we refer read- 
ers to the existing hterature (see Refs. p^-p^). 

Note that we use open boundary conditions: the blocks 
on the left and right side of Fig. |l| are not coupled to- 
gether. Although it is possible to implement the DMRG 
method with periodic boundary conditions it turns out 
that the accuracy is lower than in the open boundary 
condition case ||l5|. Therefore the method is best suited 
to study properties of two dimensional classical systems 
in contact with walls or with free surfaces. 



IV. ZERO-GRAVITY CASE 

In the limit g — > the model described by the Hamil- 
tonian (^ has been solved exactly 0| . In this section we 
summarize briefiy some known facts and present results 
obtained from DMRG calculations. 

For temperatures below a temperature Td{L), which 
Parry and Evans ^ interpreted as the critical temper- 
ature of the confined system, two phases coexist. Typ- 
ical magnetization profiles are shown in Fig. ^a). For 
Td{L) < T < Tc the system is in a single interface-like 
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state as depicted in Fig. 2(b). For L ^ oo Td{L) scales 
as 0: 

T4L)^T^^L~'/^'. (3) 

Here (3s is the exponent describing the divergence of the 
thickness of the wetting layer for a semi-infinite system 
with a surface field hi. I ^ {T^ - T)->^^. 

Swift et al. who analyzed the same system from 
a different angle, interpreted Td{L) rather as a shifted 
wetting temperature, in contrast with the point of view 
of Ref. We will come back to this point again later. 
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FIG. 2. Magnetization profiles for the model in absence of 
gravity in the two phase coexistence region for < T < Td{L) 
(a) and in the single phase region for Td{L) < T < Tc (b). 
mo denotes the bulk magnetization of the Ising model. 

We have calculated Td{L) from the correlation func- 
tion between two neighboring spins at the center of the 
strip: 

cl/2 = (sLjS^^^^j). (4) 

In the two phase coexistence region is large and 
positive since the two spins are preferably aligned. If an 
interface is present, cl/2 drops to smaller values, since in 
many configurations when the interface is located at the 
center of the strip the two spins tend to have opposite 
values. We identify Td{L) as the maximum of the tem- 
perature derivative of Cj^ ■ There is obviously no sharp 



phase transition on an L x oo system and in the present 
case Td{L) corresponds to a pseudocritical point. 

As L cxD, Td{L) scales as (|^) towards the wetting 
temperature, with the two dimensional wetting exponent 
= 1 [|13|. Figure | shows a plot of Td{L) vs l/L for 
various values of the surface field hi. On the vertical axis 
the exact values of the wetting temperatures Tw{hi), de- 
rived from Eq. (|^), are shown. As can be seen from the 
figure the scaling behavior of Td{L) is in good agreement 
with Eq. (|). 
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FIG. 3. Plot of Td{L) vs l/L for J — 1 and various values 
of the surface field hi . The crosses on the vertical axis are the 
values of the wetting temperatures derived from Eq.(^) where 
the data points are expected to scale to for L —> ca. 

The shift of Td{L) as L — > oo was also investigated in 
Monte Carlo simulations by Albano et al. [Q. They also 
found good agreement with (H). 

Table I shows a comparison between the values of the 
wetting temperatures extrapolated from a finite size scal- 
ing of the DMRG data for Td{L) and the exact values 
given by Eq. (|^); the agreement is up to the fourth dec- 
imal digit. We stress that the relevant source of errors 
for Td{L) is in the calculation of the numerical derivative 
of Ci/2; DMRG calculations provide values of 0^/2 which 
are very accurate. 



TABLE I. Comparison between DMRG data and exact values of the wetting temperatures for different values of the 
surface field h\. The number between parenthesis are the error on the last two digits of the DMRG results. 



hi 


0.1 


0.5 


0.7 


0.9 


0.99 


DMRG 


2.25761(77) 


1.95814(77) 


1.63532(77) 


1.10753(77) 


0.58848(77) 


Exact ^ 


2.25710 


1.95845 


1.63496 


1.10745 


0.58845 



V. EFFECT OF GRAVITY with all spins positive or negative are degenerate, with 

an energy per spin equal to: 



It is convenient, before discussing the full phase di- 
agram of the system described by the Hamiltonian (|^) _ 
to focus on its ground state properties. The two states ° ~ I ^ L 



1 



(5) 
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If g is large and positive the ground state configuration 
is: 



+ 1 i/ 1 < i < ■§ 
-1 i/ # + 1 < « < 



with an energy per spin equal to 



„ 3\ 2hi gL 



(6) 



(7) 



Alternatively, if g is large and negative the ground state 
has all spins reversed with respect to the configuration 
(0); in this case the energy is 



, „ 3\ ^ 2hi gL 



(8) 



The ground state is thus double degenerate if eo < 
and eo < e i-, which yields the following condition: 



4(/ti - J) 

£2 



< 9 < 



4(/ii + J) 
L2 



(9) 



At T = the two phase coexistence region shrinks as 
l/L^ towards the g = axis. The previous calculation is 
valid only for hi < J; if the surface field is larger than 
J (this implies also = 0) the range of values of g for 
which one has phase coexistence at T = is given by: 



< 5 < 



4(/ti + J) 
L2 ■ 



(10) 



In the rest of the paper we will consider only the case 
hi < J. 

The phase diagram of the model in a {g, T) plane for 
hi = 0.5, J = 1 is shown in Fig. IJ. The curves indi- 
cate the phase boundaries between the two phase coex- 
istence region (area below the curves) and a single phase 
region (above) for different values of the strip width L 
pO| . When continued down to T = the phase bound- 
aries meet the horizontal axis at the values of g equal to 
the extremes of the interval (|^). 

The phase boundaries cross the axis g ~ a,t the in- 
terface delocalization temperatures Td{L) (indicated by 
horizontal arrows in Fig. |), which scale to the wetting 
temperature T^,, indicated by a cross on the g = axis 
of Fig. ^. The striking feature of the phase diagram ||] 
is that for non-zero gravity phase coexistence extends at 
higher temperatures with respect to the case g = due to 
the competing effect of surface and gravitational fields. 
At negative g the two phase coexistence is suppressed 
at lower temperatures than for g = 0, since in this case 
gravity enhances the effect of the surface fields. 

The phase boundary maxima [(7max(-^),'7max(i)], indi- 
cated by the vertical arrows in Fig. ^, shift towards the 
bulk critical point T = Tc, g = and can be identified 
as the finite system (pseudo) critical points. As L — s- oo 
the whole critical region shifts towards the 5 = axis. 



This is in agreement with the results of Van Leeuwen and 
Sengers |11| who pointed out that in presence of gravity 
there is no criticality in an infinitely extended system. 
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FIG. 4. Phase diagram of the model in the {g, T) plane. 
The curves are the phase boundaries between the two phase 
coexistence and the single phase regions for different values 
of the strip width. Notice the scaling to the wetting temper- 
ature (horizontal arrows) and to the critical point (vertical 
arrows). 

The phase diagram of Fig. |^ is in agreement with 
the mean field results of Rogiers and Indekeu who, 
in addition, found tricritical points located in the phase 
boundaries, which separate continuous from first order 
transition lines. Obviously these features are not found 
in two dimensions where the wetting transition is always 
critical, but they should be found in three dimensions 
where the wetting tricritical point is at T > 0. 

Finite size scaling |p|~pT[| predicts that the critical point 
shifts as follows: 



T'max(^) — Tc ^ L 



Vt 



9n 



c(i) ~ L 



(11) 



(12) 



yT and yn are the thermal and magnetic exponents of 
the Ising model, which in two dimensions are yx = ^ and 
VH = 15/8. 

TABLE II. The values of the extrapolated critical temper- 
ature from the scaling analysis of Tmax(J^) for different values 
of the surface field hi, when hi < J. The exact value is 
Tc = 2.269185 . . . 



hi 


0.1 


0.5 


0.99 


DMRG 


2.269(3) 


2.272(3) 


2.271(3) 
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In deriving Eq. ( 1121) one assumes, as done by van 
Leeuwen and Sengers , that the gravitational constant 
g times a length scales as a bulk constant field. This re- 
lates the scaling of g to the bulk magnetic exponent tjh- 

The DMRG data for T^^xiL) are in very good agree- 
ment with the scaling relation (ITTl), as shown in Fig. ||. 
Table || shows the values of Tc extrapolated from the fi- 
nite size scaling analysis of T-[aax{L) by means of iterated 
fits. Results are in good agreement with the exact value. 



exponent (recall that Ai = 1/2 for the two-dimensional 
Ising model). The fact that the surface field enters in 
the form of a scaling variable hxL^^y^ is a direct con- 
sequence of the scaling form of the singular part of the 
surface free energy ]9|,p^. 

For L —* oo one should recover from ( p^ ) the scaling 
relation (|lj), which implies that: 

lim n{x) = n^^0 (14) 
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FIG. 5. Scaling of rmax(^) — c vs 1/L for different surface 
fields. To avoid overlapping of data the values of Tniax(i) 
have been shifted by a constant c, with c = for hi — 0.1, 
c = 0.3 for hi — 0.5 and c — 0.6 for hi = 0.99. Dotted lines 
are linear fits of the data points. 

The finite size scaling along the gravitational field di- 
rection is more intriguing. Figure |^ shows a plot of 
ln[5,„ax(-^)] vs ln(L) for hi = 0.5 and hi — 0.99. In 
both cases there is agreement with the scaling relation 
(|l|): as L grows the data points approach the dashed 
lines which have slope —2.875 (the value of the exponent 
— 1 — UH for the two dimensional Ising model). Notice 
that for hi = 0.99 the asymptotic behavior sets in al- 
ready for L > Lq w 20 while for hi = 0.5 this happens 
only for L > Lq ~ 60. 

The scaling behavior of (?max(^) for smaller surface 
fields is shown in Fig. 0; in this case the deviation from 
the expected exponents is so strong that one concludes 
that either ( [l^ ) does not hold for small hi or the asymp- 
totic behavior sets in only for strip widths much larger 
than those analyzed. 

Following Fisher and Nakanishi |^,|l^ who investigated 
the critical point shift in an Ising model confined between 
identical walls, one expects a scaling of gn\a.x{L) of the 
type: 



5max(i) = n {hiL^^y^) , (13) 

where fl{x) is a scaling function and Ai the surface gap 



When hi is sufficiently large and for not too small val- 
ues of L, the scaling function in ( [l3| ) can be replaced by 
its asymptotic value floo- In this limit, which we refer 
to as the saturated regime, 5max(i) becomes practically 
independent of hi . 
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ln(L) ln(L) 
FIG. 6. Plot of ln[gmax(i)l vs InL for hi = 0.5 (left) and 
hi = 0.99 (right). Error bars are smaller than symbol sizes. 
The dashed lines correspond to a slope —2.875 and are drawn 
as a guide to the eye. 

From the analysis of the numerical data for hi = 0.99 
one finds that fl{x) saturates for x > xo ~ a/20. Conse- 
quently, for hi — 0.5, 0.2 and 0.1 the saturated regime is 
expected to occur for L > Lq ~ 80, L > Lq ~ 500 and 
L > Lo ~ 2000 respectively. Notice that a saturation 
value of Lq « 80 for hi = 0.5 is in agreement with our 
numerical data. Strips of width L — 500 or L = 2000 are 
beyond the possibilities of our numerical investigation; 
actually, calculations for L > 100 (the largest size ana- 
lyzed in the present work) are feasible, but for such large 
systems the value of 5max(i) is very small and affected 
by large relative error bars that make the scaling analysis 
difhcult. 

For hi — one has 5max(-^) = since the phase bound- 
aries of Fig. 1^ become symmetric with respect to the 
g = axis. For very small surface fields one expects that 
<?max(-^) scales linearly with hi p^; this implies that: 



5 



for X — > 0. 



(15) 



Therefore, in the strongly "undersaturated" regime, 
where x = hiL^^y^ <C "\/20, one expects a scahng be- 
havior in L of the type: 



(16) 



In Fig. the dashed-dotted hnes have slopes equal to 
—2.375, which is the value of the exponent —1 — uh + 
AiUt for the two dimensional Ising model. As can be 
seen from the figure the numerical data are in good agree- 
ment with a scaling of type (16). 
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ln(L) ln(L) 
FIG. 7. As in Fig. | for hi = 0.1 (left) and hi = 0.2 
(right). Dashed and dashed-dotted lines have slopes —2.875 
and —2.375 respectively. 

We stress that ( p^ is not the asymptotic behavior of 
5max(i) as L — > oo, since for sufficiently large L the 
scaling is of type (^2|). To verify this we have calcu- 
lated 5max(i) for different strip widths at constant val- 
ues of a; = hiL^^y^ . The data points are shown in 
Fig. ^ and they are calculated in the undersaturated 
regime {x < xq ~ \/20); the agreement with the expo- 
nent I + dh = 2.875 is very good. 

To conclude this section we point out that also for the 
temperature one expects a shift of the type P,pll: 



(17) 



Nakanishi and Fisher analyzed the behavior of r(a::) 
in a mean field model confined between identical walls 
and found that r(a::) depends very weakly on its argu- 
ment X. This is also found in the present study, as it 
can be seen from the fact that (1) the slopes of the data 
points of Fig. || are almost equal, i.e. they do not de- 
pend sensibly on the value of the surface field, and (2) 
the points are very well fitted by straight lines. 
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FIG. 8. Plot of ln[gina.x{L)] vs InL at constant values of 
X = ftiL'^i^^. The data correspond to (a) x = y/TO and (b) 
X — \/0A. Dashed lines have slope —2.875. 



VI. MAGNETIZATION PROFILES 

Figures ^ and ^ show some examples of magnetiza- 
tion profiles calculated by DMRG for various values of 
the gravitational constant and temperatures. The pro- 
files (a) refer to points below the phase boundaries of 
Fig. ^, i.e. in the two phase coexistence region. The 
two coexisting phases are expected to have magnetiza- 
tion profiles similar to those depicted in Fig. ||(a); the 
profiles (a) are actually the result of the average over the 
two phases. 

The magnetization in the two phase coexistence region 
does not decay to a bulk constant value far from the walls 
as in Fig. ||(a) . The inset of Fig. || shows an enlargement 
of the profile (a) at the center of the strip: due to the 
presence of gravity the profile follows a straight slightly 
inclined line. We stress that the accuracy of the DMRG 
results for the magnetization profiles is very high, so er- 
ror bars are much smaller than symbol sizes also in the 
scale of the inset of Fig. ^. 

In the single phase region (Figs. |9pC| (b,c)) the sys- 
tem develops an interface: gravity has the tendency to 
localize the interface at the center of the strip in a con- 
figuration of minimal energy. The larger the value of \g\ 
the stronger the localization effect, as it can be seen in 
the figures: the profiles (b) correspond to a value of the 
gravitational constant 7.5 times larger than the profiles 
(c) . Notice also in the profiles (b) the competing effect of 
surface fields and gravity which is visible in the vicinity 
of the walls. 

The dashed lines of Figs. p[]lO| are the magnetization 
profiles predicted by the capillary wave theory, which are 
derived in the Appendix and are given by: 



6 



m{l) 



2mo 



dt e 



(18) 



Here 2^j_ is the average interfacial width and mp is the 
bulk magnetization of the Ising model in absence of grav- 
ity. All the parameters appearing in (1181) are known ex- 
actly, and the dashed lines of Fig. 
of a fitting. 



10 are not the results 
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FIG. 9. Magnetization profiles at a fixed temperature 
T = 2.0 and for three different values of the gravitational 
constants calculated by DMRG (filled circles). The profile 
(a) is in two phase coexistence region, (b) and (c) are profiles 
in the single phase region situated to the left and to the right 
of the two phase coexistence region of Fig. ^ respectively. 
The dashed lines are the result of the capillary wave theory 
calculation as given by (p^. Inset: enlargement of the profile 
(a) in the center of the strip. 
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FIG. 10. As in Fig. g at higher temperature (T = 2.1). No- 
tice that the discrepancy between the profiles in the interfacial 
region ((b) and (c)) and those predicted by the capillary wave 
theory increase at higher temperatures. Inset: enlargement 
of the profile (b) in the vicinity of one wall. 



Capillary wave theory profiles agree very well with the 
DMRG results especially at low temperatures where the 
approximations introduced are very good. Eq. (|l^) is 
valid in the limit ^ L where the effects of the walls 
can be neglected. This condition is satisfied at low tem- 
peratures and at not too small values of \g\. Notice 
also that if |g| is too large the magnetization profile ( p^ ) 
far from the interfacial region differs sensibly from the 
DMRG results, since the effect of gravity in that region 
has been neglected. This can be seen more clearly in the 
inset of Fig. pi. 



VII. CONCLUSIONS 



In this article we have studied the critical behavior of 
an Ising model confined between opposing walls and sub- 
ject to a bulk "gravitational" field. The competing effects 
of surface and bulk fields restore phase coexistence up to 
the bulk critical temperature, in agreement with the re- 
sults of a mean field study of the model . The strong 
thermal fluctuations in two dimensions do not affect the 
mean field results, which a fortiori should also be valid 
in three dimensions where fiuctuation effects are weaker. 

Wetting plays an important role in the model as it 
was found in the studies in the absence of gravity [0-1). 
However, limiting the analysis to 5 = causes us to miss 
much of the interesting physics that arises when gravity 
is included. In particular, one misses the critical point 
of the confined system which we have identified with 
[ffmax(-^),T'max(i)], the maximum of the phase bound- 
aries separating the two-phase coexistence from the sin- 
gle phase regions. 

We have performed a detailed analysis of the critical 
point shift as L ^ 00 and found that temperature and 
gravitational constant scale in agreement with previous 
finite size scaling hypothesis [p[-pl|. Along the gravita- 
tional field direction in the limit of small surface fields a 
crossover behavior between two different scaling regimes 
is found. 

This limit has been considered recently in studies of 

Desai et al. 



critical adsorption |22 2^. Desai et al. [g2| in a experi- 
ment on a binary liquid mixture were able to investigate 
the weak surface field regime by chemically modifying the 
surface of the solid substrate which is in contact with the 
liquid. They found an unexpected behavior of the ad- 
sorption for small surface fields, in possible disagreement 
with scaling theory. This issue was discussed recently 
also by Ciach et al. . In general we expect interesting 
physics to arise in the limit /ii — > due to the inter- 
play between bulk and surface criticality, as we found in 
the analysis of the scaling behavior of the gravitational 
constant for the model studied in this article. 
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APPENDIX: CAPILLARY- WAVE THEORY 

In capillary wave theory ||2^ one assumes that a sharp 
interface separates two regions of constant magnetiza- 
tion. We take magnetizations equal to ±mo, the bulk 
magnetization of the Ising model in absence of external 
fields, which is known exactly. In this approximation 
gravity does not affect the magnetization far from the 
interfacial region and bulk fluctuations are neglected. In 
the continuum limit the interface is described by a single 
valued function l{y) where y is the coordinate along the 
wall {—oo<y<oo) and I denotes the displacement of 
the interface from the center of the strip (see Fig, 



11) 



This is a solid-on-solid (SOS) approximation where over- 
hangs are neglected |25|. 

The continuum Hamiltonian is given by: 



H[l{y)] = 



' \ ao f dl 
''U[d-y 



U(l) 



(19) 



where ctq is the surface tension and U {I) the potential 
acting on the interface. The partition function is given 
by: 



Vl{y) e-^^['(^)] 



(20) 



where (3 denotes the inverse temperature. In (^0|) we in- 
tegrate over all the possible interface shapes described 
by single valued functions. Due to well-known relations 
between path integrals and quantum mechanics [ p6| , the 
previous problem can be mapped into a one dimensional 
quantum problem which consists in solving the following 
Schrodinger equation: 



1 



2cto/32 dp 



(21) 



The ground state wave function squared |-0o(OP denotes 
the probability of finding the interface at a position 
The potential has the following form: 



U{l;L) = Wil;L) + Vg{l) 



(22) 



Here W{l;L) is the confining potential due to the pres- 
ence of the walls: 



W{l;L) 



if|/|<L/2 
-foo if|/|>L/2 



(23) 



Vg{l) represents the contribution of gravity to the interfa- 
cial potential and can be calculated from the microscopic 
Hamiltonian ([T]). If the interface is located between the 
k-th and (k+l)-th spin of the system gravity gives the 
following contribution to the energy: 



Q(fc) = gmo I ^(2i - 1 - L) - ^ (2i - 1 - L) 



i=k+l 



2mogk{k — L) 



(24) 



Shifting appropriately the origin of the coordinates one 
finds: 



Vg{l) = 2mogf + c 
where c is a constant. 



(25) 




FIG. 11. Example of an interfacial configuration de- 
scribed by a continuous single- valued function l{y) as in the 
solid-on-solid approximation. All the spins at the two sides 
of the interface are fixed and equal to ±mo. The arrows de- 
note the elastic force that tends to bring the interface in its 
equilibrium position. 

Neglecting the effect of the confining potential ( p3| ) be- 
comes the Schrodinger equation for a harmonic oscillator, 
with ground state wave function equal to: 



Ml) 



denotes the interfacial width: 



6 



T 



(26) 



(27) 



Here toq is the bulk magnetization and do is the sur- 
face tension which are known exactly for the d = 2 Ising 
model a.t g — 0. The analysis is valid for ^± <C i where 
the probability of finding the particle outside the walls 
is negligible and the confining potential W{l;L) can be 
ignored. 

From (^6|) one can easily calculate the magnetization 
profile as function of the distance from the center of the 
strip I: 
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771(0 ^rno{ I ds \Ms)\^ - r ds |i^o(s)|' I (28) 



which yields the result given in (nl 
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